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Nucleon Decay Matrix Elements with Nf = and 2 Domain- Wall Quarks 

Yasumichi Aoki^ [RBC collaboration] * 

'^Physics Department, University of Wuppertal, Gaussstr 20, 42119 Wuppertal, Germany. 

The nucleon decay matrix elements of three-quark operators are calculated with domain-wall fermions. Opera- 
tors are renormalized non-perturbatively to match the MS (NDR) scheme at NLO. Quenched simulation studies 
involve both direct measurement of the matrix elements and the chiral Lagrangian parameters, a and (3. We also 
report on the dynamical quark effects on these parameters. 
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1. INTRODUCTION 

Nucleon decay once observed, is the hall- 
mark of the physics beyond the standard model. 
On-going deep mine experiments, though yet to 
observe an event, are pushing up the lower bound 
of the lifetime of the nucleon [2|, excluding GUT 
models which allow nucleons decay more fre- 
quently. 

The dominant decay mode of the nucleon is to 
a pseudoscalar meson and a lepton, where low 
energy hadronic interactions are important. The 
factorization technique, commonly used for the 
hadron decays, leads to the decay amplitude writ- 
ten in terms of the Wilson coefficients and the 
low-energy matrix elements of the dimension-six 
operator consisting of three quarks and one lep- 
ton. Various QCD model calculations had esti- 
mated the hadronic part of the matrix element 
with results varying over an order of magnitude 
(see a compilation in ^). The amplitude must be 
squared in the width or lifetime, meaning two or- 
ders of magnitude different lifetime from different 
estimations. The initial lattice calculations show 
smaller but, still large deviation (see a compila- 
tion in Hj). 

Recent lattice calculations concentrate on how 
to remove the systematic uncertainties, especially 
the discretization error, by either taking the con- 
tinuum limit 5 of the Wilson fcrmion or em- 
ploying the domain- wall fermion (DWF) |HI7| . 
All these calculations, however, have been done 
within the quenched approximation. 
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We report here the updated results of the nu- 
cleon decay matrix elements with DWF with 
quenched approximation as well as the estimate of 
the low energy parameters with Nf = 2 dynam- 
ical DWF. Detailed calculations will be given in 
the forthcoming publication 0. 

2. QUENCHED SIMULATION 

The parameters of our quenched simulation are 
/3 = 0.87 for the gauge coupling of the DBW2 
gauge action, domain- wall height M5 = 1.8, fifth 
dimension Lg — 12, the lattice volume 16^ x 32 
where latter is the temporal direction. The in- 
verse lattice spacing ~ 1.3 GeV is obtained 
by the nip = 0.77 GeV input at the chiral limit, 
which gives the physical spatial lattice size as 
L ~ 2.4 fm. While we only have one lattice spac- 
ing, the scaling violation is expected to be small 
for DWF as the other hadronic observables have 
shown 9 . 

The renormalization of the operator is done in 
combination of the non-perturbative renormaliza- 
tion (NPR) with the continuum perturbation the- 
ory m NLO. The NPR scheme UHI is employed for 
the lattice operator to match the MOM scheme. 
The one-loop matching factor that we have cal- 
culated converts it to the MS in NDR scheme. 
Then, using the two-loop anomalous dimension 
|11| . the operators are run to /i = 2 GeV. 

We demonstrated the NPR work well last year 
reusing the data for the quark bilinear renor- 
malization at Lg = 16. We have performed NPR 
calculation at Ls ~ 12, where the matrix elements 
are obtained as well, to get the proper renormal- 
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Figure 1. Summary of Wq (Eq.^ with direct and 
indirect method in quenched calculation. Opera- 
tors are renormalized at fi = 2 GeV. 

ization factor and to confirm the basic properties 
such as the absence of the mixing. 51 configura- 
tions are used for NPR. 

We also showed the lattice value of the hadronic 
matrix elements in the previous reports |6I7| . 
While the number of independent gauge configu- 
rations analyzed in this study is unchanged (100 
configurations), we make use of discrete symme- 
try transformation properties which relate the dif- 
ferent Green functions to be averaged. In this way 
the effective statistics have been increased . Fur- 
ther we have changed the analyses details to get 
more robust results. All these have changed the 
results slightly and helped shrinking the error. 

In Fig. n we show the results of the relevant 
form factor Wq, which, for example in the p — > 
71^^/+ decay, defined through 



7r;p\0%/^.Jp; k) = Pl[Wo - %m^W 



(1) 



where of/^.^ = e'3Hu'^CPR/Ld^)PLu\ is 
the momentum transfer from proton (fc^) to pion 
(PAt)j 111 = — p^, Up is the proton spinor. The 
results both with direct and indirect calculations 
are shown. The former involves various two- and 
three-point functions with many different param- 
eter values to allow extrapolation to the physical 
kinematics point, thus is expensive. The latter is 
obtained with help of tree-level chiral perturba- 
tion theory |12I4| where the low energy parame- 



ters a and f3, defined as 

aPLUp = (0|<Jp), (3PlUp = (Olojjp), (2) 

are calculated on the lattice. The other low en- 
ergy parameters are taken from the experiment 
[?]• The indirect method uses only a few two 
point functions, thus less computational effort is 
required than the direct method. The SU(2) fia- 
vor symmetry of u and d quarks relates the dif- 
ferent matrix elements, e.g., 



-plL:L\P)=^iA0l,.Ap). 



(3) 



Every other possible matrix element is identical 
to one of those in Eq.^or Fig.Qup to sign factor. 

The final tt state matrix elements has apparent 
deviation between the direct and indirect calcu- 
lations. This may simply mean the limit of the 
tree-level xPT, as the pion has large momentum. 
Indeed in the soft pion limit (p^ — s- 0) the ex- 
pected relation (WQ — iqiWq)^'''" = olI\/2]t^ holds 
numerically within the error. We note that the fi- 
nal K state matrix elements are estimated better 
than those for the tt by the indirect method. 

The values of individual matrix elements are 



different from those obtained for 



2.3 GeV 



with Wilson fermion 0], while ratios of the ma- 
trix elements are similar. We consider our DWF 
results are closer to the continuum limit. Indeed 
it is the case for the low energy parameters. 

The results of a and /3 are summarized in the 
Table n] They have a different relative sign and 
the relation a -I- /3 = approximately holds with 
1.3(7. Recent quenched calculation results 5._ with 
Wilson fermion in the continuum limit are also 
shown. There the operators are renormalized by 
an improved lattice perturbation theory. We re- 
mark that the DWF results at coarser lattice are 
consistent with them. 

3. DYNAMICAL QUARK EFFECTS 

To investigate the quenching error we calculate 
the a and /? parameters on the TV/ = 2 dynamical 
DWF configurations ^Hj . The simulation param- 
eters are the same as quenched, except for gauge 
couphng /3 = 0.8, which gives ~ 1.7 GeV, 
L ~ 1.9 fm by mp in the chiral limit. We have 
three dynamical quark masses {rridyn) spanning 
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Figure 2. Pion mass dependence of a in both 
dynamical [N j — 2) and quenched (N j — 0) cal- 
culations. 

Tfisll ^ TTidyn ^ "T-s, where is the strange 
quark mass. The spatial lattice size divided by 
the pion Compton wave length is Lm^ — 4.7 at 
the smallest. We expect the systematic error from 
the finite volume effects is subdominant given the 
large statistical error of the matrix elements. The 
operator renormalization is done in the same way 
as quenched calculation. The number of configu- 
rations used is ~ 100 for matrix elements or ~ 40 
for NPR for each dynamical mass. 

The measurements are done at the valence 
masses equal to the dynamical masses. The pion 
mass dependence of |a| for both quenched and dy- 
namical calculations are shown in Fig.|2| The dy- 
namical result has stronger dependence than 
quenched. After rather long extrapolation to the 
chiral limit with linear function of quark mass we 
obtain the a and /3 parameters as shown in the 
Table [T| The relation a + (3 = holds within 
the error. The values are consistent with the 
quenched ones. We note that these estimates 
from lattice QCD lie in the middle between the 
smallest (0.003 GeV^ [HI) and largest (0.03 GeV^ 
estimates of QCD model calculations. 

4. OUTLOOK 

The present dynamical DWF simulation with 
quark masses rudyn ~ has not revealed the 

unqucnching effect on the a and (3 parameters 
in the chiral limit with 20% statistical error. To 



Table 1 

Summary of the low energy parameter of nucleon 
decay at the renormalization scale fi — 2 GeV. 
Quoted errors for DWF are statistical only. 



Fermion 


Wilson 5 


DWF 


Nf 
a [fm] 
|a| [GeV3] 
|/?| [GeV3] 





0.0090(09)(li9) 


2 
0.15 0.12 
0.010(1) 0.012(2) 
0.011(1) 0.012(2) 



reduce the error to 10% level, much lighter dy- 
namical quark region needs to be explored. As to 
this direction, the direct calculation of the matrix 
elements becomes more important. Furthermore, 
test of scaling to the continuum limit will also be 
needed. 
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